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Ww\rick/  '  ABSTRACT 

Static  stability  is  investigated  as  a  contribution  to  numerical 
baroclinic-zone  analysis,  developed  by  Renard  and  Clarke,  and  presently 
being  produced  on  a  hemispheric  basis  by  the  U.S.  Navy  at  the  Fleet 
Numerical  Weather  Facility,  Monterey,  California  (FNWF) .    Significant  - 
level  and  mandatory-level  data  from  the  North  American  region  at  00Z  and 
1200Z,  20  August  1965,  are  utilized  in  an  experiment  relating  850,  700, 
500,  and  300-mb  static  stabilities  and  associated  parameters  to  a  FNWF 
front-location  parameter  proportional  to  horizontal  shear  of  the  thermal 
wind  which  is  dependent  on  mandatory  level  data  only. 

Three  methods  are  used  to  compute  static  stability  in  the  form  of 
Mcl©   ,  the  most  promising  one  involving  a  parabolic  fit  of  potential 
temperature  to  the  four  mandatory  levels .    A  derivative  of  static  stability 
is  used  to  locate  baroclinic  zone  boundaries.    Finally,  the  vertical  gradient 
of  potential  temperature  in  combination  with  the  slope  of  isentropic  surfaces 
is  used  to  determine  isobaric  gradients  of  potential  temperature  which 
yields  an  analysis  of  baroclinicity  with  resolution  surpassing  that  of  the 
present  FNWF  product. 
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1 .        Introduction 

Considerable  developmental  work  has  been  performed  at  the  U.S. 
Navy  Fleet  Numerical  Weather  Facility  (FNWF)  and  the  U .  S.  Naval  Post- 
graduate School  (NPGS)  on  the  problem  of  locating  frontal  zones  by 
numerical  methods  .  /  1 ,  8 ,  9_/  The  most  promising  and  currently  operational 
approach  utilizes  the  second  derivative  of  potential  temperature  (-^r&z) 
as  a  front-location  parameter,  hereafter  referred  to  as  GG0.    Axes  of  relative 
minima  and  maxima  in  the  GG0  fields  which  bound  the  baroclinic  zones 
are  associated  with  relative  maxima  in  the  horizontal  shear  of  the  thermal 
wind.    The  axes  of  relative  maxima  represent  the  warm-air  boundary  of 
baroclinic  zones  and,  as  such,  are  a  numerical  form  of  conventional 
fronts .    These  numerical  fronts  have  a  close  resemblance  to  hand-analyzed 
fronts ,  as  may  be  seen  in  Figure  1  which  shows  a  typical  FNWF  operational 

analysis  of  GG9 ,  and  Figure  2  which  indicates  the  relation  of  the  numerical 

1 
front  (i.e.  ,  GG9  ridge)  to  the  sea-level  pressure  analysis. 

However,  recent  statistical  evaluations  indicate  a  need  for. improving 
the  numerical  product  in  certain  areas,  such  as  in  eastern  ocean  and  sub- 
tropical regions  and  frontal  occlusions.    The  author  of  this  study  compared 


Throughout  the  subsequent  discussion  the  following  terminology 
is  used  interchangeably:    (a)  GG0  fronts  ,  numerical  fronts ,  baroclinic 
zone  boundaries;  (b)  GG0  ridge,  numerical  front,  warm  boundary  of 
baroclinic  zone;  (c)  GG9  trough,  cold  boundary  of  baroclinic  zone; 

(d)  numerical  frontal  zones  ,  baroclinic  or  hyperbaroclinic  zones; 

(e)  hand-analyzed  fronts,  manually  analyzed  fronts,  conventional 
fronts,  NMC  fronts. 


11 


1000-mb  fronts,  specified  by  GG0 ,  with  hand-analyzed  surface  fronts  of 
the  National  Meteorological  Center,   (NMC)  for  16  observation  times  in  the 
late  summer  of  1965  / 4/.    It  was  found  that  72%  of  the  GG9  fronts  could 
be  associated  with  NMC  fronts.    The  average  distance  separating  the  two 
fronts  was  about  90  nautical  miles,  slightly  greater  in  sparse  data  areas. 
These  statistics  were  based  on  a  sample  of  over  250,000  miles  of  fronts. 
Table  I  recapitulates  this  and  other  findings. 

In  its  present  operational  form,  0  is  determined  from  mandatory  level 
data  only,  however  0  or  its  derivative  (G0  and  GG0)  may  be  computed  along 
any  level  of  meteorological  interest  (e.g.,   1000-mb,  earth's  surface,  5,000 
feet,  mean  level  of  1000/700-mb  layer.)    It  is  reasonable  to  expect  that 
consideration  of  the  detailed  temperature  structure  (i.e.  ,  from  significant 
as  well  as  mandatory  levels)  in  the  vertical  will  contribute  to  the  accuracy 
of  locating  baroclinic  zones  at  the  aforementioned  levels  . 

This  study  introduces  detailed  vertical  resolution  through  a  static 
stability  parameter  computed  from  significant-level  as  well  as  mandatory- 
level  data.    The  quantity  was  evaluated,  first  as  an  entity  unto  itself,  and 
later,  combined  with  other  quantities  to  determine  frontal  locations  and 
intensities.    Specifically,  the  computed  static  stability  field  was  considered: 

(1)  Relative  to  FNWF's  operational  model  of  GG0,  as  well  as  to  their 
present  static  stability  evaluation  /5 ,  6/.     (See  section  4) 

(2)  As  a  contribution  toward  a  more  detailed  GG0  analysis  through 

the  introduction  of  detailed  vertical  temperature  structure,     (see  section  6) 

2.         Selection  of  a  Suitable  Static  Stability  Parameter. 

Static  stability  refers  to  the  resistive  force  offered  to  motions  of 
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various  types  by  the  vertical  density  structure  of  the  atmosphere  in  the 
gravitational  field  and,  specifically,  for  the  atmosphere  in  hydrostatic 
equilibrium.    Here  static  stability  is  proportional  to  o&  -  X  .    When    &.>  i    , 
an  individual  air  parcel  displaced  dry  adiabatically  in  a  vertical  direction 
will  tend  to  return  to  its  original  position,  and  the  air  column  is  said  to  be 
statically  stable. 

Frontal  layers  are  zones  of  relatively  high  values  of  static  stability, 
since  within  such  layers  temperature  usually  deviates  from  its  normal 
decrease  with  increasing  height  to  yield  isothermal  or  inversion  conditions. 

Thus ,  it  appears  meaningful  to  compute  static  stability  for  geometric 
or  isobaric  levels  of  interest  over  an  area.    Relatively  high  values  of  static 
stability  ought  to  be  located  within  frontal  zones  and,  thus,  the  static 
stability  should  have  a  high  positive  correlation  with  baroclinicity;  GG9 
ridges  and  troughs  should  bound  zones  of  high  static  stability.    Since  a 
large  value  of  stability  does  not  imply  the  existence  of  a  baroclinic  zone, 
and  vice-versa,  it  is  expected  that  some  stability  ridges  (that  is,  axes  of 
relative  maxima)  and  numerical  fronts  may  exist  independently.    Figure  3 
shows  a  low  tropospheric  vertical  cross-section  of  the  eastern  Pacific  ■/_!_/. 
This  cross-section  illustrates  high  values  of  static  stability  and  baro- 
clinicity in  the  east  Pacific  trade-wind  inversion  east  of  135W,  while  also 
indicating  an  area  of  poor  correlation  between  baroclinicity  and  stability  in 
the  inversion  west  of  135W. 

Gates  has  made  an  analysis  of  advantages  and  disadvantages  of 

several  forms  of  static  stability  /8/.    Selection  of  an  optimum  stability 

parameter  for  this  project  was  based  on  several  self-imposed  requirements: 
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(1)  static  stability  will  be  computed  on  constant  pressure  surfaces   ; 

(2)  the  measure  should  have  equal,  or  nearly  equal,  values  at  different 
levels  when  %     is  the  same  fraction  of     od      ;  (3)  static  stability  should 
be  zero  when    X   =  ^i     ;  (4)  input  parameters  should  be  readily  available 
from  radiosonde  data;  and,  (5)  the  computation  must  readily  lend  itself  to 
both  manual  and  numerical  analysis,  the  latter  currently  existent  or  planned 
for  the  future .    Requirements  1,4,  and  5  are  self  explanatory .    Requirements 
2  and  3  are  necessitated  by  a  desire  to  have  comparative  fields  of  static 
stability  at  the  levels  investigated. 

A  relatively  simple,  non-dimensional  measure  possessing  the  desirable 
properties  is: 

->^  ~      e       dp  M 

This  measure  meets  the  above  requirements  very  well,  moreover,  it -Is 
directly  proportional  to     oA.   -  &    .    Gates  found  that  SS  computed  by  this 
method  showed  a  slight  decrease  with- height  up  to  300-mb,  followed  by  an 
increase,  but  the  range  of  this  variation  is  small  and  will  not  be  of  any 
consequence  here.    Equation  (1)  also  portrays  the  stratosphere  realistically, 
and  shows  the  relatively  unstable  upper  troposphere  in  both  summer  and 
winter. 
3.        Area  and  Time  of  Experiment. 

The  area  and  time  selected  for  computing  SS  were  all  of  the  United 
States,  Canada,  and  Alaska  at  00Z  and    122  \  20  August  1965.    The  reported 

2 

The  present  FNWF  operational  frontal  model  utilizes  mean  9  and 

its  derivatives  for  the  layer  1000/700-mb  as  the  basic  frontal  parameter. 
Derivatives  are  computed  at  constant  pressure.    (See  Figures  1,  2) 

14 


radiosonde  temperature  data  were  then  utilized  in  three  different  methods 
to  compute  values  of  SS ,  for  each  of  approximately  100  reporting  stations, 
at  850,  700,  500,  and  300-mb. 

20  August  1965  was  a  well  documented  day,  with  extensive  data 
readily  available  from  listings  provided  by  FNWF.    The  date  was  marked 
by  the  presence  of  a  strong  cold  front  and  large  stability  over  the  eastern 
and  southern  United  States.    As  Gates  pointed  out,  wintertime  situations 
give  greater  stability  contrasts  (up  to  50%) ,  especially  in  the  lower  tropo- 
sphere,  and  better  delineate  frontal  zones,  but  20  August  proved  to  be  quite 
satisfactory  for  a  pilot  project. 
4.         Computation  and  Areal  Analysis  of  Static  Stability. 

I.      Methods  for  Computing  Static  Stability. 

a.    Method  A.    This  method  consists  of  computing  the  stability 

3 
for  the  two  significant  layers     adjacent  to  a  mandatory  level  and  plotting 

the  larger  of  these  two  values  on  the  appropriate  constant  pressure  chart 

(850,  700,  500,  and  300-mbs.)    This  approach  has  the  desirable  effect  of 

sharply  defining  the  stable  baroclinic  zone  boundaries  by  assigning  to  the 

mandatory  level  the  highest  possible  SS  value  from  one  of  the  two  adjacent 

significant  layers.    Figure  4  shows  a  portion  of  a  typical  sounding  for  which 

SS  might  be  computed  using  Method  A.    Since  layer  AB  is  more  stable  than 

layer  BC ,  SS  is  computed  for  AB  and  this  value  is  assigned  to  500-mb,  the 

mandatory  level.    Here,  the  500-mb  level  represents  the  warm-air  boundary 

A  significant  layer  is  bounded  by  two  significant  levels,  one  or 
both  of  which  may  be  mandatory  levels  as  well. 
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of  a  stable  and,  presumably,  barociinic  zone.    Thus  a  maximum  stability 
value  found  in  a  layer  within  the  stable  zone  is  placed  at  its  boundary. 
Stability  within  the  significant  layer  was  computed  according 
to  the  definition  in  Equation  (1) .    P  and  9  are  the  median  values  for  the 
layer,  while  the  layer  value  of     i£     is  approximated  by    £§     of  the  layer. 
In  this  method,  a  significant  layer  was  defined  as  being  bounded  by  a 
mandatory  level  and  its  nearest  significant  level,  regardless  of  the 
layer's  depth.    In  another  method  to  be  discussed  later,  layers  of  less 

than  10-mb  thickness  were  ignored. 

4 
A  Fortran  60  program    was  written  to  compute  stability  to  four 

significant  figures  by  this  method.    The  values  were  rounded  to  three 

significant  figures  and  plotted  and  contoured  on  FNWF's  1:30  million  polar 

steoreographic  projection.    Resulting  analyses  appear  as  Figures  5-8,  with 

comments  on  the  technique  appearing  in  Appendix  I . 

Figure  6  (see  station  206)  displays  the  merit  of  this  approach, 
since  large  gradients  of  stability  do  allow  the  stable  zone  boundaries  to 
be  easily  located  in  this  case. 

But  Method  A  was  not  considered  satisfactory  for  several  reasons 
First,  it  must  be  remembered  that  the  objective  is  to  compute  stability  at 
the  mandatory  levels.    This  system  arbitrarily  assigns  a  layer  value  as  the 
mandatory-level  value.    The  actual  value  at  the  mandatory  level  approaches 
this  arbitrary  value  as  a  limit,  but  for  many  stations  this  method  gave  the 

4This  and  all  subsequent  programs  were  run  on  the  NPGS's 
Control  Data  Corporation  1604  digital  computer. 
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stability  for  a  small  layer  near  the  mandatory  level,  which  is  not  nec- 
essarily representative  of  stability _at  the  mandatory  level. 

Secondly,  this  method  introduced  vertical  variation  into  hori- 
zontal gradients  of  SS .    Figure  9A  shows  a  hypothetical  case  in  which  SS 
is  higher  in  the  significant  layer  above  station  A  than  in  the  significant 
layer  below  station  B.    The  resulting  field  of  SS  tends  to  give  excessive, 
unrealistic  gradients  of  SS  within  the  baroclinic  zone,  as  well  as  suggesting 
improper  boundaries  of  the  zone.    This  problem  is  magnified  for  stable  zones 
which  are  thin  (less  than  20-mb  in  depth)  and  undergo  marked  variations  in 
depth  over  small  horizontal  distances. 

Thirdly,  it  will  be  noted  from  a  discussion  of  Methods  B  and  C, 
that  the  latter  give  a  more  meteorologically  satisfying  field  of  SS  relative 
to  GG0  baroclinic  zones  and  the  hand-analyzed  fronts  of  NMC. 

Consequently  this  method  did  not  look  as  promising  as  Methods 
B  and  C,  and  was  not  further  developed. 

b.    Method  B  also  utilized  Equation  (1)  for  each  of  the  four 
mandatory  levels  of  interest  by  employing  a  mathematical  best  fit  of  9  at 
the  mandatory  level,  using  transmitted  radiosonde  data.    This  data  consists 
of  temperature  and  pressure  recorded  at  each  significant  level  (level  where 
a  significant  change  in  temperature  lapse  rate  occurs) ,  as  well  as  for  all 
mandatory  levels.    Thus  the  sounding,  as  plotted  from  transmitted  data, 
is  discontinuous  at  each  of  these  levels.    However,  analysis  of  radiosonde 
data,  assuming  every  temperature  contact  is  a  significant  level,  indicates 
the  transmitted  discontinuity  points  are  largely  manufactured  by  coding  and 
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processing  procedures.    Thus,  fitting  a  continuous  curve  to  the  transmitted 
sounding  for  the  purpose  of  obtaining  **   is  reasonable  and  justified  /_2j . 
This  is  accomplished  by  numerically  solving  the  following  system  of 
simultaneous  linear  equations  for  each  sounding: 


Ou.  -    a  +  bpu-+  cp*.z 

0vnL=   a -i-  bpmL+cpmi.2  '2) 

9l  -  a+  bfu  -t-cpL2' 


A  Fortran  60  program  using  a  Jordan  2  subroutine  was  written  to 
solve  this  system  for  the  constants  a,  b,  and  c,  since  all  six  pressures 
and  temperatures  (hence,  0)  are  known  from  the  sounding. 

At  the  mandatory  level,  the  vertical  gradient  of  potential  tem- 
perature, -^—  becomes 

5©rnL       _       k  +  ?C.JDrf\U  (3) 


±-  -     b  +  ZC-P^l 


Thus  ,  Equation  (l)  was  numerically  evaluated  for 


SSmL--^    x-|^b      -_-|^(b+zcpmL)  (4) 


It  can  be  shown  (Appendix  III)  that  the  finite  difference  form 
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of     4^     is: 

7^y$t*^  +i^*Ae*]  (5) 

This  expression  was  used  to  manually  check  computerized 
results  of  Equation  (4) .    One  such  computation  appears  in  Appendix  II. 


For  FNWF  numerical  operational  purposes,  Equation  (5)  would  undoubtedly 
prove  to  be  a  more  desirable  system  than  Equation  (3) . 

Computed  static  stabilities  were  plotted  and  contoured  at  a  .05 
interval,  similar  to  values  resulting  from  Method  A. 

Two  problems  appeared  at  this  point:    (a)  the  existence  of  some 
negative  values  of  SS  on  each  chart  (i.e.  ,  absolutely  unstable  conditions) 
which  verified  to  be  real  and  representative  of  low  level  summer  conditions, 
and, thus,  were  accepted;  and  (b)  a  few  questionable  values  which  caused 
irregularities  in  otherwise  smooth  contour  patterns.    Re-examination  of 
sounding  data  revealed  that  a  majority  of  these  questionable  values  of  SS 
resulted  when  a  significant  level  was  located  within  10  mb  of  a  mandatory 
level . 

When  thin  layers  of  a  few  millibars  are  considered,  normal 
errors  in  pressure  and  potential  temperature  can  vary  the  value  of  SS  at 
the  mandatory  level  by  a  large  percentage.    This  problem  is  discussed  more 
fully  under  Method  C. 

Furthermore,  layers  of  a  few  millibars'  thickness  may  be  con- 
sidered mesoscale  in  character  and  not  sufficiently  deep  to  define  synoptic 
scale  frontal  layers,  or,  alternately,  as  nearly  horizontal,  thin,  stable 
layers  whose  relation  to  baroclinic  zones  is  unexpectedly  poor.    For 
example,  near  700-mb,  a  10-mb  layer  (about  400  feet)  representing  a 
narrow  baroclinic  zone  100  miles  wide,  yields  a  1/1500  slope.    Since 
baroclinic  zones  as  depicted  by  GG0  average  almost  300  miles,  a  thin 
layer  (10  mb  or  less)  is,  for  all  practical  purposes,  horizontal.   (Figure  10) 
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By  ignoring  layers  of  10  mb  or  less  ,  this  system  would  be  throwing  out 
only  the  nearly  horizontal  stable  layers . 

A  third  method  has  now  been  suggested,  one  which  does  not 
consider  layers  of  less  than  10-mb  thickness  adjacent  to  a  mandatory  level. 
This  is  Method  C,  discussed  below.    Since  Method  A  altered  58  out  of  a 
possible  440  values  computed  in  Method  B,  taking  all  four  levels  into 
account,  the  location  of  SS  ridge  lines  differed  for  each  method,  but  contour 
patterns  were  basically  the  same.    For  this  reason,  the  contour  analysis  by 
Method  B  was  not  included  in  this  paper.    Instead,  Figures  11-14  compare 
the  locations  of  the  SS  ridge  lines  produced  by  Methods  B  and  C.    These 
figures  are  discussed  in  Appendix  II.    On  each  chart,  SS  values  are  shown 
for  all  stations  where  SS  varied  depending  upon  the  method  used.    The 
number  plotted  above  the  station  is  the  value  computed  by  Method  B;  the 
number  below  the  station  is  Method  C's  result.    It  is  evident  that  the  choice 
of  methods  heavily  influences  the  location  of  ridge  lines,  and  the  choice  of 
a  minimum  significant-layer  thickness  can  be  a  critical  factor  in  locating 
SS  zones . 

c.    Method  C .    This  method  is  considered  as  best  representing 
SS  at  the  four  levels  of  interest.    It  is  identical  to  Method  B,  with  one 
exception.    Layers  10  mb  and  less  in  thickness  are  considered  unrepresent- 
ative.   In  cases  where  a  significant  level  exists  within  10  mb  of  a  mandatory 
level,  the  next  lower  or  higher  significant  level  is  substituted.    For  example, 
assume  the  levels  of  interest  are  886,  858,  850,  and  795  mb.    For  850-mb 
analysis,  858  mb  would  be  ignored  since  it  is  within  10  mb .    Instead  of 
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858  mb,  886  mb  and  its  associated  temperature  are  substituted  for  P 

and  T    . 
u 

Thus,  Method  B  is  a  very  precise  method,  taking  every  fluc- 
tuation of  temperature  in  the  layers  surrounding  the  mandatory  levels  into 
account.    Method  C  includes  a  small    "smoothing  factor",  which  has  the 
effect  of  filtering  out  some  of  the  "noise"  created  by  very  thin  (i.e.  ,  non- 
synoptic  scale)  layers.    Furthermore,  the  minute  temperature  changes  in 
thin  layers  dictates  use  of  layers  thicker  than  10  mb ,  primarily  because  of 
the  accuracy  demanded  by  this  method. 

When  compared  with  FNWF  numerical  fronts  for  the  same  time 
and  level,  the  more  satisfying  appearance  of  the  SS  ridges  using  Method  C 
raises  the  question  as  to  the  best  smoothing  factor  to  be  employed.    Per- 
haps raising  the  minimum  thickness  of  the  significant  layer  to  20  mb  may 
even  be  justified.    An  error  analysis  of  SS  brings  out  this  point.    Consider- 
ing an  error  in  static  stability  (ASS)  at  a  mandatory  level  as  a  function  of 
errors  in  significant  layer  depth  /  A  (A  P]_/,  and  temperature  change  in  the 
layer  /  A  (A  ©)_//  Equation  (1)  may  be  differentiated  to  obtain: 

ASSmt  "-  -L-(A9)fA#  A(Af)+VAfAS)]^  (6) 

Probable  errors  in  pressure  and  temperature  are  assumed  to  be  1  mb  and 
1°C,  respectively,  in  the  low  troposphere  /l_l/ '.    So  as  not  to  inflate  the 
proposed  error,  it  may  be  assumed  the  probable  error  occurs  at  only  one 
boundary  of  the  significant  layer,  errors  in  0  and  T  are  the  same  magnitude, 
and  transmission  errors  are  excluded.    Further,  assuming  P  =  700-mb, 
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9  =  310°,  and     A0  =  .  2°C  for  each   AP  interval  of  2  mb,  the  following 
results: 

AP(mb)  &e(°C)  &SS 

2                                  .2  104% 

10                                    1  21% 

20                                  2  11% 

Thus,  dramatic  changes  in  accuracy  are  achieved  by  proper 
selection  of  the  smoothing  factor,  while  still  retaining  the  desired  vertical 
resolution.    Further  experimentation  is  certainly  warranted  before  final 
selection  of  the  optimum  minimal  thickness  of  significant  layers. 

d.    Analysis  Techniques  Employed  with  Method  C  Data . 
Data  obtained  by  Method  C  were  plotted  on  1:30  million 
polar  stereographic  charts  and  analyzed  for  00Z  and  T2Z      20  August  at  all 
four  levels.     (Figures  15-24) 

f 

Any  attempt  to  contour  an  area  as  large  as  the  United  States , 
Canada,  and  Alaska  on  the  basis  of  about  100  observations  is  open  to 
criticism,  since  any  one  of  several  interpretations  may  be  used  without 
violating  data.    This  is  especially  true  in  Canada,  where  extensive  sparse 
data  areas  exist.    The  proper  orientation  of  the  isolines  of  SS  is  intrinsic 
to  proper  location  of  the  ridge  line  of  SS .    Since  the  ridge  line  suggests  the 
location  of  the  frontal  zone  and/or  its  boundaries,  an  assist  or  first  guess 
was  sought  for  the  SS  analysis. 

Since  the  ultimate  goal  of  this  project  is  to  guide,  or  even 
replace,  the  present  FNWF  system  of  frontal  analysis,  it  was  considered 
valid  to  use  the  GG0  field  as  an  aid  in  contouring  SS  by  orienting  the  SS 
lines  parallel  to  the  numerically  computed  baroclinic  zone  boundaries 
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wherever  data  allowed.    It  turned  out  that  forcing  was  not  necessary, 
since  SS  isolines  did  indeed  parallel  baroclinic  boundaries  on  most 
occasions,  especially  in  areas  of  strong  baroclinicity .    The  result  looks 
very  promising,  since  the  patterns  obtained  show  a  striking  resemblance 
to  GG0  ridge  lines,  notably  at  850  and  700  mb.     (Figures  15,  16,  20,  and 
21) 

II.     Discussion  and  Evaluation  of  SS  Charts  Produced 

At  this  juncture  it  was  felt  that  Method  C,  considering  baroclinic 
zones  presented  SS  in  its  most  palatable  form  as  representative  of  frontal 
phenomena.    The  only  obvious  shortcoming  is  its  subjectivity,  being  a 
manually  produced  analysis  product. 

Appendix  III  contains  a  detailed  analysis  of  eight  charts,  one 
for  each  level  of  interest  for  both  00Z  and  1200Z  20  August  1965,  the  charts 
appearing  as  Figures  15-18  and  20-23.      Speaking  in  general  terms,  it  can 
be  observed  that: 

a.  Most  ridge  lines  of  SS  are  in  close  proximity  to  NMC  fronts 
and  to  ridge  lines  of  GG0;  exceptions  are  manifest  at  higher  levels  (500  and 
300-mb) ,  and  in  mountainous  areas. 

b.  There  are  more  miles  of  SS  ridge  lines  than  GG9  fronts,  and 
more  GG0  fronts  than  NMC  fronts,  the  latter  fact  indicated  in  /_!_/  and  /4/. 

c.  SS  ridge  lines  have  obvious  vertical  consistency  between  the 
850  and  700-mb  levels,  especially  over  the  dense  data,  non-mountainous 
region  of  North  America,  where  a  hand-analyzed  front  also  exists.    Here, 
the  700-mb  front  trails  its  850-mb  counterpart  by  about  100  miles.    This 
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seems  consistent  with  the  normal  slope  of  a  low  tropospheric  front.    In 
the  same  area,  the  500-mb  SS  ridge  trails  the  700-mb  position  by  as  much 
as  300  miles  (Great  Lakes  area,   1200Z)  and  it  is  at  this  level  that  vertical 
continuity  breaks  down.    Not  all  fronts  reach  the  500-mb  level,  and  of  those 
that  do,  many  show  a  sharply  decreased  slope.    Such  an  observation  is  in 
agreement  with  findings  by  other  investigators  /10_/.    Vertical  consistency 
of  SS  ridges  in  mountainous  regions,  sparse  data  areas,  and  the  peripheral 
map  area  is  fair  to  poor. 

d.  The  ridge  lines  at  300-mb  are  questionable  representatives 
of  low  troposphere  fronts  for  the  same  reasons  outlined  above.    There  is  the 
additional  problem  posed  by  the  tropopause .    In  computing  SS  by  Method  C , 
data  up  to  250-mb  were  often  used.    In  winter,  the  low-level  tropopause 
gives  high  stabilities  over  some  mid-latitude  stations,  and  thus  gives  a 
false  indication  of  a  front.    Several  soundings  in  Canada  (stations  722, 
909,  913,  915,  917,  and  several  others)  indicated  a  tropopause  existing 
between  250-"and JOO-rnb 

e.  The  ridge  lines  of  SS  showed  good  time  continuity  in  the  area 
of  good  vertical  consistency. 

f.  In  areas  of  sparse  data,  subjectivity  enters  the  analysis.    In 
areas  where  three  or  four  reports  determine  the  contour  pattern  over  an  area 
of  thousands  of  square  miles,  baroclinic  zones,  history,  and  vertical  con- 
sistency are  the  only  aids  available.    These  are  important  considerations 
when  viewing  application  on  a  hemispheric  scale,  utilizing  computerized 
solutions  . 
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III.     Comparison  of  SS  with  FNWF's  Static  Stability  Analysis  ( (T  ) . 

a .  Observations: 

FNWF's  data  processing  program  produces  mean  static  stability 
for  five  layers  (0"  )  /5,6/.    Such  a  computation,  of  course,  was  not  meant 
for  aid  in  frontal  analysis.    However,  it  is  pertinent  to  this  investigation  to 
compare  the  merits  of  SS  computed  for  a  level,  using  the  maximum  vertical 
resolution,  with  CT,  computed  for  a  layer  using  mandatory  level  information 
only.    A  detailed  comparison  of  the  two  analyses  appears  in  the  comment 
section  of  Appendix  IV.    The  discussion  here  is  limited  to  the  following 
general  observations. 

The  use  of  mandatory  level  data  only  gives  a  gross  picture  of 
stability,  and  says  little  about  stability  at  a  particular  level.    A  frontal 
layer  could  be  embedded  within  the  775/600-mb  layer,  and  mandatory  level 
data  might  never  detect  it.    On  the  other  hand,  SS  using  significant  as  well 
mandatory  level  data  is  more  representative,  if  for  no  other  reason  than  the 
increase  in  data  and  the  vertical  resolution. 

SS  has  the  added  advantage  of  being  computed  for  four  specific 
constant  pressure  surfaces  in  the  atmosphere.    With  the  wide  use  of  850}  , 
700f  ,  and  500-mb  charts,  stability  computed  at  a  level  is  more  meaningful 
than  stability  within  a  gross  layer. 

b.  Conclusions: 

When  compared  with  FNWF's  C  analysis,  SS  is  a  more 
meaningful  representation  of  frontal  zones  because  of: 
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(1)  the  use  of  actual  data  from  radiosonde  reports  vice 
processed  data  as  used  by  FNWF. 

(2)  the  high  vertical  resolution  resulting  from  the  use  of 
significant  vice  mandatory  level  data. 

(3)  the  improved  accuracy  derived  from  using  a  parabolic  fit 
of  potential  temperature  at  the  mandatory  level  with  no  assumption  about 
the  vertical  gradient  of  static  stability. 

IV.    Recommendations  on  SS  . 

Based  on  the  aforementioned  findings,  the  following  recommenda- 
tions are  made  relative  to  further  testing  and  application  to  frontal  analysis: 

a.  Limit  the  mandatory  levels  to  the  850,    70CT>  ,  and  500-mb 
surfaces  for  low  tropospheric  frontal  phenomena. 

b.  Compute  SS  as  indicated  by  Method  C. 

c.  Statistically  compare  SS  ridges  with  computer  and  hand- 
drawn  fronts  and  baroclinic  zones  produced  by  NMC  and  FNWF,  respectively 

d.  Consider  programming  computation  of  SS  on  a  hemispheric 
basis,  say  in  the  framework  of  the  data  processing  scheme  at  FNWF.    It  is 
realized  that  this  requires,  apriori ,  a  suitable  data-handling  program  in 
order  to  marry  the  mandatory  and  significant  level  radiosonde  data.    Pre- 
sently, the  significant  level  data  are  copied  on  tape  but  not  read  out  or 
utilized  in  any  way  in  FNWF's  automatic  data  processing  program. 

e.  Subsequent  to  d,  analyze  the  SS  field  considering  it  or  its 
derivatives  as  an  adjunct  to  the  GG9  analysis  for  locating  baroclinic  zones. 
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f .    Investigate  the  feasibility  of  computing  SS  for  the 
surface,  or  1000  mb,  by  assigning  the  SS  computed  for  the  lowest 
significant  layer  (or  level)  as  the  value  at  the  surface. 
5.        The  Gradient  of  Stability:    GSS  . 
I.      Definition  of  Parameter. 

As  indicated  in  the  Introduction,  FNWF's  present  operational 
product  for  frontal  analysis  is  the  GG9  field.     This  field,  as  defined  by 
Renard  and  Clarke,  is  the  directional  derivative  of  the  gradient  of  potential 
temperature  along  its  gradient  /_&_/.    In  mathematical  terms: 

GG0  =  -V/v*0/-$§r=-V/7©/-r£©  5  ~-$r*  (7) 


Referring  to  SS  and  recalling  that  it  is  already  a  first  derivative, 
a  similar  expression  involving  stability  may  be  written: 

GGS=  -VSS-^fj-     E    -VSS  .7?*  (8) 

The  GG0  and  GSS  fields  are  entirely  similar  in  gross  features, as 
may  be  seen  by  Figure  25. 

II.     Discussion  of  Procedure. 

GSS  was  calculated  manually  for  the  700-mb  level  for  00Z 
20  August  1965.    The  mechanics  of  this  calculation  were  as  follows: 

a.  Using  FNWF's  hemispheric  printout  of  potential  temperature 
at  700-mb,  rf.    was  determined  at  each  grid  point  over  the  area  of  interest. 

b.  For  each  grid  point,  the  gradient  of  stability  was  deter- 
mined using  -Tnfer-    /  where     6fn_  was  taken  as  120  nautical  miles,  with 

A  n©  © 

the  grid  point  at  the  centrum.    This  distance  is  rather  arbitrary,  being 
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selected  mostly  on  the  basis  of  allowing  for  a  well  defined  field  of  GSS. 
In  other  words,  the  larger  the  distance  An0/  the  flatter  the  field  of  GSS. 

c.    GSS  was  plotted  on  a  1:30  million  polar  stereographic 
projection  and  then  analyzed.    This  chart  appears  as  Figure  25.    The  warm 
air  boundary  of  a  baroclinic  zone  is  associated  with  the  ridge  or  axis  of  a 
relative  maximum  in  GSS,  while  the  cold  air  boundary  is  representative  of 
the  trough  or  axis  of  a  minimum  is  GSS . 
III.    Conclusions. 

a.  GSS  does  appear  to  offer  some  advantage  over  the  SS  field 
for  frontal  information.    SS  defines  a  ridge  only,  which  is  associated  with 
the  centrum  of  the  zone,  while  GSS  defines  a  ridge  and  a  trough  which  out- 
line the  boundaries  of  the  baroclinic  zone.    In  comparing  these  boundaries 
with  the  baroclinic  zone  limits  derived  from  GG0,  the  GSS  boundaries  define 
a  zone  which  averages  about  100  miles  narrower  than  GG9,  except  over  the 
west  central  U.S.     (See  Figure  25)    There  are  at  least  two  reasons  for  the 
suggested  narrower  baroclinic  zones: 

(1)  greater  vertical  resolution  in  computed    GSS  vs.  GGQ 
from  FNWF;  or 

(2)  use  of  120  miles  to  compute  GSS  instead  of  several  grid 
lengths,  as  used  by  FNWF. 

b.  The  manual  computation  of  GSS  is  a  time  consuming  process, 

requiring  about  three  hours  to  produce  the  700-mb  chart  for  the  North 

American  area.    However,  the  computation  of  GSS  is  adaptable  to  the 

computer,  the  mechanics  of  the  operation  being  almost  identical  to  the 

computation  of  GG0. 
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IV .    Recommendations . 

a.  If  GSS  is  investigated  at  some  later  date,  the  procedures 
outlined  in  this  section  are  recommended,  but  it  is  suggested  that  the 
entire  operation  be  programmed  with  GSS  as  the  final  product. 

b.  If  the  above  recommendation  is  carried  out,  statistically 
compare  GSS  with  GG0  to  determine  if  the  former  is  consistently  better 
in  outlining  the  baroclinic  zones. 

c.  The  following  section  concerns  an  alternate  approach  to 
obtaining  baroclinic  zone  boundaries.    Since  SS  and  GSS  experiments  were 
but  pilot  projects ,  it  is  recommended  that  a  more  exhaustive  study  be  made 
relating  the  GSS  scheme  to  the  approach  in  the  next  section. 

6.         G0  Computed  Using  High  Vertical  Resolution. 
I.      Definition  of  Parameter. 

As  discussed  in  Section  3,  GG0  in  its  present  operational  form 
utilizes  mandatory  level  data  only  for  frontal  analysis. 

Consider  the  following: 

If      Q~-4c\n.  (f,T\&) 

d-G^i^+lf^p  (9) 

Next,  set  d0;O  for  purposes  of  computing  the  slope  of  an 
isentropic  surface,  here,  in  the  direction  of  maximum  slope.    Thus, 


i£ 


$n< 


50      /P6 


•-    fW#  "                                          do) 

»*•-£.=  Hfrl*-3r  (11) 

where,    GO  *  |— 8fU  <12) 

Finally,  GG9  =    -G(G0)  =  -  ^y^z 
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II.     Discussion  of  Procedure  . 

Equation  (11)  can  be  solved  from  available  data  since  %^p 
may  be  computed  using  Method  C,  and    |-|£"l      can  be  determined 
utilizing  FNWF's  hemispheric  printout  of  P    .    Thus,  a  G9  results  which 
incorporates  the  high  vertical  resolution  inherent  in  the  use  of  significant 
level  data  through  consideration  of    ^    .    This  new  G9  may  then  be  com- 
pared with  FNWF's  hemispheric  printout  of  G0  to  determine  whether  high 
vertical  resolution  used  in  this  fashion  enhances  the  operational  frontal 
product. 

The  mechanics  of  determining  Igj^L  involved  obtaining  hemi- 
spheric printouts  of    P    for  each  2 . 5  A  of  potential  temperature  from  280A 
to  330A,  in  order  to  compute   -3£>    at  850 ■'-  and  700-mb.  for  00Z  20  August 

on© 

1965.     |  ^f)Q|©    was  computed  at  each  station,  using  a  centered  finite 
difference  scheme  applied  to  two  mesh  lengths  of  the  FNWF  grid.    Inter- 
polation of    |-§^|Q  for  the  exact  0  surface  was  made,  as  appropriate. 


Finally,  the  products  of    \4£-      and  g%     were  obtained  for  each  station 

at  85u>    and  700-mb.    The  resulting  values  of  GO  (hereafter  referred  to 

as  computed  GO)   were  transferred  to  1:30  million  charts  and  compared  with 

FNWF's  GO  values. 

III.    Results. 

The  comparison  between  computed  and  FNWF  GO  values  appears 

in  tabular  form  as  Appendix  V.    Note  that  computed  GO  has  been  used  as  the 

basis,  since  its  values  are  correct  to  three  significant  figures,  and  its  use 

allows  the  selection  of  177  specific  reference  points  (considering  both  charts) 

for  comparison  of  the  values  . 
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Each  chart  was  arbitrarily  separated  into  two  parts,  for  reasons 

which  will  soon  become  evident.    One  part  was  the  baroclinic  zone,  which 

o  5 

was  defined  here  as  the  area  enclosed  by  the  1.0  C/100  km.  isoiine.      The 

remaining  area  was  simply  considered  outside  the  baroclinic  zone.    Figures 

26  and  27  show  a  comparison  of  the  baroclinic  zones  from  computed  and 

FNWF  G0  values. 

The  table  in  Appendix  V  reveals  that  computed  G0  gives  relatively 
higher  values  in  baroclinic  zones  and  lower  values  outside  the  zones.    This 
fact  is  obvious  when  one  notes  that  of  51  reference  points  within  the  baro- 
clinic zone,  G0  is  less  than  -.2  C/lOOkm  at  two  points  only;  while  of  126 
reference  points  outside  the  baroclinic  zone,  only  three  values  of  t\  G0 
are  above  +  .2°C/l00km. 

Thus,  as  shown  by  the  medians  obtained  for  each  area,  computed 
G0  gives  a  wider  range  of  values  overall,  implying  a  sharper  change  in 
gradient  (i.e.  ,  GG0)  across  the  boundaries  of  the  zones. 

The  16  observations  at  850-mb  in  which  the  FNWF  G0  exceeded 
computed  G0  by  more  than  1.0  C/lOOkm  fell  mostly  in  the  southwest  desert 
area  of  the  United  States,  as  may  be  inferred  from  Figure  26.    As  stated  in 
Section  1,  this  area  poses  a  problem  for  numerical  frontal  analysis  by  GG0 
in  that  fronts  are  indicated  where  none  exist  in  manual  analyses,  especially 
in  the  summer.    Computed  G0    has  markedly  low  values  in  this  area.    Thus, 
either  the  vertical  resolution  of   ^    or  the  use  of  a  three-point  vice  a 

Values  on  both  charts  ranged  from  nearly  0  to  3.0  C/100  km. 
with  two  isolated  higher  values  on  computed  G0:  5.15°C/l00km  at  station 
764  and  3.50°C/100  km  at  station  518. 


31 


five-point  derivative  in  computing     -|£    has,  at  least  for  this  map  time, 
resolved  an  area  of  difficulty. 

A  final  result  concerning  locations  of  relative  maxima  and 
minima  in  the  field  show  that  of  23  relative  maxima  on  the  computed  G0 
field,  there  were  17  maxima  on  FNWF's  printout  which  fell  within  one  mesh 
length.    Of  29  relative  minima,  FNWF's  printout  indicated  23  minima  within 
one  mesh  length.    These  facts,  combined  with  a  cursory  glance  at  Figures 
26  and  27,  show  a  similarity  of  contour  patterns,  except  over  the  southwest 
U.  S.  at  850-mb. 

IV.    Conclusions. 

It  is  difficult  to  make  any  sweeping  statements  about  G0 
computed  by  Equation  (11)  on  the  basis  of  one  map  time,  but  the  following 
comments  certainly  appear  justified: 

a.  There  is  a  resemblance  in  the  shape  and  location  of 
baroclinic  zones  defined  by  the  two  methods  ,  especially  the  warm  air 
boundary . 

b.  Computed  G0  defines  a  narrower  baroclinic  zone  with 
higher  values  inside  and  lower  values  outside  the  zone. 

c.  Computed  G0  seems  to  eliminate  the  troublesome  area 
over  the  mountainous  and  desert  southwest  at  850-  and  700-mb. 

Thus,  computed  G0  appears  to  show  considerable  promise, 
especially  in  light  of  b.  and  c.    The  question  reduces  to  an  analysis  of 
Equation  (11)  and  a  decision  as  to  whether  the  vertical  resolution  obtained 
by  the  use  of  significant  level  data  has  narrowed  and  sharpened  the 
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baroclinic  zone,  or  whether  the  use  of  a  three-point  vice  a  five-point  deri- 
vative is  responsible.    Comparisons  made  by  Associate  Professor  R.J. 
Renard  of  the  U.S.  Naval  Postgraduate  School  and  Leo  C.  Clarke  of 
FNWF,  Monterey,  California,  revealed  that  the  five-point  method  actually 
gives  a  stronger  gradient  than  the  three-point.    Thus,  it  certainly  appears 
that    -^     plays  a  dominant  role  in  statement  b . 

As  a  test,  three  stations  at  700-mb  were  arbitrarily  selected 
and  a  comparison  made  between  the  three-point  (two  mesh  lengths) 
derivative  used  in  this  experiment  and  the  five-point  technique  presently 
used  operationally  by  FNWF.    The  five-point  derivative  is  described  by  the 
quartic  interpolation  polynomial 

-fft9AT  i*o*  r£a  £Ww  +  a£crt!*«»^>fm2-><W'*i3  (13) 

Where  d=  mesh  length  =  381km.  at  60     latitude,  i  =  0  is  the  station 
location,  i  +  1  and  i  +  2  are  locations  one  and  two  mesh  lengths  from  the 
station,  respectively,  in  a  direction  normal  to  isobars  on  a  constant  0 
surface.    The  table  below  suggests  that  the  gradients  computed  by  the  two 
methods  are  similar,  but  certainly  not  the  same. 

Sp/  8  ne(mb/l00km) 
Station  three  -points  five-points 

712  32.8  29.3 

722  40.6  45.9 

311  22.3  14.6 

V .      Recommendations . 

Continue  experimentation  on  computed  G0  by: 
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a .  programming  the  computer  to  produce   ~%r\&    from  PG 
charts  and  multiplying  the  result  by  i^t    from  Method  C; 

b.  operate  on  computed  G0  to  get  GGO  (this  is  already  an 
operational  procedure  at  FNWF)  with  high  vertical  resolution; 

c.  make  a  detailed  statistical  comparison  of  the  two  GG0 
parameters  to  determine  which  product  best  delineates  the  baroclinic  zone 
and  its  associated  frontal  phenomena. 

7.         General  Conclusions  and  Recommendations. 

In  view  of  the  stated  conclusions  in  each  section,  only  a  brief  overall 
review  will  be  given  here. 

This  entire  study  utilizes  data  from  two  map  times  and  attempts  to 
show  the  value  of  using  significant  as  well  as  mandatory  level  data  as  a 
contribution  toward  numerical  baroclinic  zone  and  frontal  analysis.    Three 
different  parameters  for  frontal  analysis  have  been  presented:    SS ,  GSS ,  and 
a  computed  G9.    It  is  not  the  purpose  of  this  investigation  to  select  from 
these  three  a  single  proposal  which  is  considered  best.    All  three  have 
strong  and  weak  points ,  and  a  statistical  study  covering  several  periods  of 
time  in  at  least  two  seasons  would  be  required  to  make  an  intelligent  selec- 
tion . 

It  does  appear,  however,  that  the  detailed  thermal  structure  of  the 
atmosphere  in  the  vicinity  of  the  mandatory  levels  can  be  better  visualized 
by  use  of  significant  level  data,  and  that  the  computation  of    ^^      using  a 
parabolic  fit  is  a  most  promising  approach. 
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TABLE  1 

EVALUATION  ANALYSIS 

FNWF  (NUMERICAL)  1000-MB  vs  NMC  (MANUAL)  SFC  FRONTS 

Dates:    16  obsn.  times  in  period  23  Aug  -  6  Sept  65 

FNWF  FRONTS:    259,000  n. mi. 

NMC  FRONTS:    284,000  n. mi. 

%  FNWF  FRONTS  assoc.  with  NMC  FRONTS:        72% 

%  NMC  FRONTS  assoc.  with  FNWF  FRONTS:        66% 


Sector 

Frontal 
Distance 
(n .  mi . ) 

Average 

(FNWF 

Summer 

Difference  (°lat) 
vs  NMC) 

Winter 

North  America 

49,250 

1.6 

1.4 

Atlantic 

46,900 

1.4 

1.8 

Pacific 

46,300 

2.1 

2.3 

Europe-North  Africa     26,000 
Asia  18,000 


1.8 
1.2 


See/1/ 
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APPENDIX  I 
SS  Analysis  by  Method  A 
Technique  and  Criteria. 

SS  computations  by  Method  A  were  plotted  on  1:30  million  polar 
stereographic  paper  and  the  resulting  contour  analysis  appears  as  Figures 
5-8.    Values  of  SS  vary  from  slightly  negative  to  +.6,  with  most  values 
falling  in  the  +.05  to  +.3  range.    Contour  interval  has  been  arbitrarily  set 
at  .05  for  values  from  0  to  .2,  and  at  .1  thereafter. 

Contours  of  SS  are  generally  elongated  and  elliptical,  somewhat 
resembling  a  jet  stream  analysis.    The  elongation  generally  appears 
parallel  to  ridge  lines. 

Ridge  lines  are  drawn  only  for  relative  maximum  values  greater  than 
.  1  and  then  only  if  more  than  one  station  shows  such  a  magnitude .    Ridges 
are  also  omitted  for  areas  of  weak  gradient  (although  values  of  SS  may 
exceed  .1)  as  well  as  when  considerable  subjectivity  is  required  in  locating 
the  ridge.    Although  this  discussion  concerns  Method  A,  the  principles 
apply  to  every  analysis  attempted  in  this  experiment. 
850-mb  00Z  20  August  1965.     (Figure  5) 

It  will  be  recalled  from  Section  4  that  Method  A  has  the  advantage 
of  producing  numerical  values  of  SS  which  are  somewhat  higher  than  those 
computed  by  the  other  methods  discussed  herein.    These  higher  values  have 
produced  a  ridge  over  the  eastern  U.S.  which  is  well  defined  by  high  SS  at 
stations  712  y  528,  520,  429,  and  403.    This  ridge  has  not  been  connected 
to  the  ridge  over  the  central  U.S.  because  of  low  SS  at  station  327,  but  a 
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slight  relaxing  of  criteria  could  have  allowed  such  a  marriage.    The  other 
three  ridges  on  this  chart  are  a  function  of  a  very  few  station  reports  ,  and 
exact  placement  was  difficult. 

The  western  U.S.  has  few  contours  because  many  of  these  stations 
are  at  pressures  greater  than  850-mb.    The  zero  contours  in  the  U.S. 
southwest  desert  region  bounds  areas  of  absolutely  unstable  lapse  rates. 

The  overall  pattern  has  some  resemblance  to  the  frontal  patterns 
achieved  by  other  methods  (see  Figure  15),  but  it  is  clear  that  correlation 
is  too  poor  for  a  comparative  analysis  to  be  of  any  value. 
700-mb  00Z  20  August  1965.     (Figure  6) 

This  chart  is  characterized  by  several  high  values  of  SS  which  dictate 
the  location  of  ridge  lines .    Most  of  these  values  are  computed  from  thin 
layers  adjacent  to  the  700-mb  level,  and  the  overall  pattern  shows  very 
poor  correlation  with  patterns  produced  by  later  methods  and  operational 
products,     (see  Figure  16)    The  most  striking  example  appears  at  station  206 
Other  computations  modify  this  1.86  reading  to  a  value  of  about  .1,  a 
remarkable  change.    The  same  criticism  can  be  made  of  stations  722,  907, 
and  913. 

With  the  exception  of  the  ridge  over  the  eastern  U.S.,  vertical 
continuity  with  the  850-  and  500-mb  levels  is  poor. 
500-mb  00Z  20  August  1965.     (Figure  7) 

The  500-mb  chart,  when  compared  to  GG9  and  later  methods  in  this 
study  (Figure  17),  has  the  highest  correlation  of  any  chart  of  the  series, 
showing  a  long  ridge  of  SS  from  SE  Canada,  across  the  central  U.S., 


39 


then  north  into  central  and  northern  Canada.    High  SS  values  at  stations 
520,  712,  747,  and  913  dictate  the  position  of  the  ridge,  just  as  they  do 
in  Figure  15.    The  ridge  over  the  eastern  seaboard  is  well  defined,  but 
its  position  is  not  verified  by  any  other  method.    A  large  portion  of  the 
western  U.S.  exhibits  moderately  high  values  of  SS ,  but  the  area  is  so 
broad  and  the  gradient  so  flat  that  it  is  not  possible  to  draw  a  realistic 
ridge  line . 
300-mb  00Z  20  August  1965.     (Figure  8) 

The  pattern  at  300-mb  is  random  and  irregular.    The  only  ridge 
showing  vertical  continuity  with  lower  levels  is  the  ridge  over  eastern 
Canada  and  the  east-central  U.S.    As  indicated  in  the  text,  there  is 
evidence  to  support  the  claim  that  the  300-mb  analysis  may  not  be  an  aid 
to  lower  troposphere  frontal  analysis.    Its  value  as  a  front-locator  is, 
therefore,  questionable. 
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APPENDIX  II 
Comparison  of  SS  Ridge  Lines , 
Method  B  vs  Method  C 
Technique  and  Criteria. 

The  general  discussion  of  Methods  B  and  C  points  out  that  the 
techniques  are  the  same  except  that  in  Method  B  all  significant  layers  were 
considered. 

Figures  11-14  point  out  the  differences  in  location  of  SS  ridges  caused 
by  differences  in  SS  values  at  individual  stations  where  thin  layers  (<10mb) 
were  present.    Each  analysis  was  accomplished  without  referring  to  charts 
produced  by  the  alternate  method,  and  therefore  are  as  objective  as  possible. 

Criteria  for  ridge  location  are  the  same  as  those  discussed  in  Appendix 
I  for  Method  A. 
850-mb  00Z  20  August  1965.     (Figure  11) 

This  chart  illustrates  how  changes  in  two  values  (stations  327  and 
451)  can  adjust  the  location  of  SS  ridge  lines  by  as  much  as  200  miles. 
The  lower  value  of  .000  calculated  by  Method  B  forces  the  ridge  well  north 
of  station  327,  and  away  from  the  NMC  front  through  the  area  (Figure  15) . 
The  other  ridges  are  identical,  since  no  other  thin  layers  were  encountered. 
700-mb  QQZ  20  August  1965.     (Figure  12) 

Numerous  thin  layers  at  700-mb  force  a  Method  B  pattern  which  differs 
considerably  from  Method  C  and  GG0  ridges  as  shown  by  Figure  16.    As  can 
be  seen  in  that  figure,  Method  C  ridges  have  very  high  correlation  with  the 
FNWF  numerical  fronts,  but  Method  B  produces  a  wavy,  random  pattern 
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which  does  not  correlate  well  with  GG0.    It  appears  that  the  elimination 
of  thin  layers  (Method  C)  has  smoothed  and  adjusted  Method  B  into  a 
meaningful  SS  ridge  analysis.    Certainly  this  700-mb  chart  seems  to 
justify  the  arguments  presented  in  section  4. 
500-mb  00Z  20  August  1965.     (Figure  13) 

The  500-mb  analysis  shows,  once  again,  how  changes  in  SS  values 
tend  to  alter  the  SS  ridge  lines.    Since  verification  of  both  Methods  B  and 
C  with  FNWF's  numerical  fronts  is  fair  to  poor  on  this  chart  (see  also 
Figure  17) ,  comments  will  be  limited  to  pointing  out  that  Method  C  again 
seems  to  smooth  out  the  ridge  pattern,  yielding  a  more  meteorologically 
satisfying  pattern. 
300-mb  00Z  20  August  1965.     (Figure  14) 

Method  C  simply  modifies  a  "no  correlation"  pattern  to  another 
"no  correlation"  pattern  over  central  Canada.    Only  five  values  differ,  and 
the  pattern  by  both  methods  is  subjective  due  to  sparse  data  and  a  wide 
variation  of  values . 
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Derivation  of  Equation  (5) 


IP:      ©=a  +  bp  +  cp2 
DO     .    ,  -  _ 

r~=  D  +  2CP 

Op 

p 

©^  =  a  +  bp-^  +  cp^ 
©2  -a  +  t>P2  +  °P2 


0o=a  +  DP3  +  cp-5c 
0i-.02  =  b(p1-p2)+o(p12-p22) 

°2-°3  =  b(p2-p3)+C(P22-p32)  f(pt-p^)(p1+p2)] 

(0l-G2)p2-p3  =  b(pl-p2)(p2-p3)+C^Pl2^P22)(p2-p3? 

(o2-e3)p1.p2  =  b(p1.p2)(p2-p3)+cb22-p32)(pi.p2)] 

'"(©l  -  ©2^pu  -(e2~G9^Pl  =  c  P^PiXPi  +^)-"pu  "  ?^pu<p2  +  p3^pli] 

(©x  -  ©2^D22  "  p32)  =bW  -  P?)(^22  -p32> 
(©2  -  ©3)(px2  -  P22)  =b(P2  "  "°j)^i2  "  D22) 

(©! -©2^pu(p2 +  t>3)  "  C(Q2  "  Q3^T)l(pl +p2^    =  " 

bf$?l&Pu ( p2  +  P3))  - £& ou  ap1  ( px  +  p2 )3 J 

RECALL:      ||  =  b + 2cp 

0G       &QifcPv,(p?  +P?)  -fo&toh  (Pi+p^A  .^jAQ-,Apr  -A©  &p,3  rj 

3U  -APfS^fTn*?*^^  P£rl  fcPu^l4-^  "?AiAPl<P2  +  P3#J 

APuAPi^P  -APiAPuAP 

=  A9-.  (d,  +  t)Q  -  AQ..(p«  +  p,Q    -   2pAQi    +  2pK9„ 
APXAP  APu&p         APAPi         APAPU 

*APU  =   (P2"p3^    ^Pls=(Pl   -  P2^»    Same   for  e« 


A2 


SINCE  p  =  p2  for  mandatory  levels: 

/D0\  =  _1_  ^  /^4-p  A9  -2-o2AQ)   J_  /gp^  -  i^AQ  -  P^OX 

I^pA  ap  I =     JAPi      '  *p  v      /.^u     / 

THUS:    /30\    -    JL  /(pT-p^)^  +    (pg  -Pi?AOu J 

Up/ 2"  &p  V  ^  Api        ^pu   / 


As  an  example  of  this  equality,  consider  station  520 

Pittsburgh,  Pennsylvania,  at  850-mb. 

PRESSURE  P. TEMPERATURE 

883  298.6 

850  300.1 

771  301.8 

Computer  result,  using  parabolic  fit,  SS  =  .1082 
By  finite  differencing: 

OP   AP  L  ^pl      ^pu  J 

=  _L.  ^19  x  1.5   33  x  l.?1  =  .0383 
112  L   33        29   J 

SS  =  p.  30  =  850   x  .0383  =  .1083 
0  ^p   300.1 

Error  in  fourth  place. 
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APPENDIX  III 
SS  Analysis  by  Method  C 
Comparison  of  SS  with  GG0  and  NMC  Fronts 
Technique  and  Criteria. 

SS  computations  by  Method  C  were  plotted  and  analyzed  applying  the 
principles  discussed  in  SS  Analysis  by  Method  A.     (Appendix  I) 
850-mb  00Z  20  August  1965.     (Figure  15) 

Four  distinct  SS  ridge  lines  appear:    (1)  from  just  north  of  Newfoundland 
through  New  England  south  to  West  Virginia  and  west  to  Kansas;  (2)  from 
North  Dakota  across  western  Canada  to  the  Pacific;  (3)  extreme  northern 
Canada;  and  (4)  western  Alaska .    The  position  of  (1)  seems  well  defined  by 
a  dense  data  network  in  its  northeastern  portion  where  the  position  is 
dictated  by  reports  at  only  five  stations:  816,  815,  807,  712,  and  600. 
Since  these  five  stations  cover  250,000  square  miles,  location  northeast 
of  station  712  is  questionable.    The  position  at  712  looks  very  good,  since 
a  strong  maximum  occurs  there  at  700-mb,  indicating  that  the  700-mb  front 
is  over  the  station  and,  therefore,  the  8  50-mb  front  should  be  southeast  of 
712,  as  plotted.    Along  the  eastern  seaboard  the  ridge  line  lies  within  the 
.150  contour,  and  its  western  end  is  anchored  by  high  readings  at  stations 
456  and  451.    This  ridge  verifies  well  with  the  NMC  front,  with  differences 
of  less  than  fifty  miles  along  50%  of  its  route.    The  position  of  the  ridge 
line  seems  to  split  the  difference  between  the  NMC  front  and  the  GG9  ridge. 
The  ridge  over  western  Canada  is  based  on  just  a  few  readings,  and  would 
not  have  been  drawn  in  except  for  the  reading  at  station  867.    It  is   in- 
teresting to  note  that  the  ridge  line  could  have  been  drawn  from  station  945 
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north  to  the  strong  maximum  over  station  938.    This  points  out  the  subjec- 
tivity inherent  in  hand  analysis  in  sparse  data  areas.    The  ridge  in  northern 
Canada  is  based  entirely  on  the  enormous  maximum  over  station  938.    The 
ridge  is  well  defined  in  the  heart  of  a  baroclinic  zone.    Lack  of  data 
prevented  continuing  this  ridge  farther  east.    Negative  stability  layers  are 
found  at  several  western  stations.    These  negative  values  result  from  one 
significant  layer  showing     tf  >  ^    ,  a  not  uncommon  occurrence  over  the 
desert  areas  in  summer  at  850-mb.    It  is  also  interesting  to  note  high  SS 
along  the  California  coast.    This  is  probably  the  trade- wind  inversion 
(Figure  3)  affecting  stability  at  850-mb.    Finally,  many  readings  in  the 
Rocky  Mountain  area  are  missing  because  of  altitude. 
700-mb  00Z  20  August  1965.     (Figure  16) 

The  pattern  here  shows  four  ridges.    The  ridge  across  the  eastern 
United  States  is  the  same  ridge  which  appeared  on  the  850-mb  chart  and 
shows  excellent  vertical  consistency,  in  that  it  trails  the  850-mb  ridge 
by  about  100  miles.    If  one  considers  the  slope  of  a  cold  front  to  be  about 
1/100  between  850    and  700-mb,  the  ridge  position  looks  reasonable. 
SS  trails  GG9  over  the  eastern  one  third  of  its  path  and  the  distance  between 
the  two  ridges  broadens  to  150  miles  over  Missouri  and  Iowa,  but  northwest 
of  South  Dakota  the  two  ridges  are  in  very  close  proximity.    The  ridge  in 
northeastern  Canada  is  based  on  one  large  reading  at  station  090.    This 
ridge  could  just  as  well  have  been  omitted,  since  there  are  no  hints  as  to 
orientation  and  direction.    A  ridge  over  the  Pacific  Northwest  is  well  defined, 
but  does  not  lie  in  a  baroclinic  zone,  and  does  not  have  a  counterpart  at 
the  850-  and  500-mb  levels  .    A  ridge  across  Alaska  east  and  south  into 
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central  Canada  compares  fairly  well  with  a  GG9  ridge.    The  SS  ridge  is 
much  longer  on  its  southern  end,  however,  based  on  high  SS  over  stations 
945  and  848.    The  vertical  consistency  of  this  ridge  is  not  outstanding, 
but  this  could  well  be  a  function  of  the  sparse  data  area.    The  western 
portion  of  this  frontal  structure  which  lies  in  a  baroclinic  zone,  shows 
vertical  consistency,  but  over  central  Canada  there  is  no  baroclinic  zone, 
and  the  front  meanders  noticeably.    A  general  feature  of  this  chart  is  the 
flat  gradient  over  the  southeast,  south  central,  and  Rocky  Mountain  states. 
SS  is  not  high  except  over  station  304,  and  since  SS  increases  to  the  north, 
no  ridge  was  drawn  even  though  the  500-mb  analysis  shows  a  ridge  in  the 
vicinity. 
500-mb  00Z  20  August  1965.     (Figure  17) 

The  ridge  pattern  at  500-mb  shows  four  ridges.    The  500-mb  ridge 
over  the  eastern  United  States  shows  high  positive  correlation  with  GG9 
and  with  the  SS  ridge  at  850V    and  700-mb  from  Newfoundland  to  Ohio. 
Northwest  of  Ohio  SS  moves  steadily  away  from  GG0 ,  taking  a  more  northerly 
route  out  of  the  baroclinic  zone.    It  may  well  be  that  the  SS  ridge  does  follow 
GG0  and  maintains  vertical  continuity  with  the  ridges  at  lower  levels ,  but 
high  readings  over  stations  747  and  645,  plus  low  readings  at  945  and  879, 
make  it  difficult  to  draw  a  ridge  across  southwestern  Canada.    A  well  defined 
ridge  over  the  southern  states  is  not  considered  frontal,  since  it  is  not 
vertically  consistent  with  a  lower  level  front  and  does  not  lie  in  a  baroclinic 
zone .    The  same  comments  can  be  made  about  a  ridge  from  Colorado  to 
Washington.    However,  the  ridge  over  Alaska  into  Canada  is  vertically 
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consistent  with  ridges  mentioned  at  lower  levels ,  and  lies  in  a  baroclinic 
zone.    Whether  this  ridge  actually  ties  into  the  ridge  over  the  midwestern 
United  States  is  conjecture.    It  is  not  surprising  that  certain  difficulties 
are  encountered  at  500-mb.    Some  cold  fronts  may  not  extend  to  this 
altitude ,  or  the  character  of  the  front  may  change .    As  pointed  out  by 
Gates  ,  the  value  of  static  stability  may  begin  to  reduce  at  this  level,  a 
fact  which  may  affect  the  pattern . 
300-mb  00Z  20  August  1965.     (Figure  18) 

This  chart  appears  to  be  of  little  value  in  location  of  ridge  lines 
which  represent  fronts.    It  does  seem  that  SS  values  are  lower  at  this  level, 
as  evidenced  by  the  extensive  flat  gradient  over  the  entire  United  States, 
with  average  values  of  .060.    There  are  three  definable  ridge  lines  on  the 
chart,  but  none  are  associated  with  lower  level  ridges  and  only  the  one 
over  the  desert  southwest  lies  in  a  baroclinic  zone. 

It  seems  logical  that  fronts,  if  they  exist  at  all,  will  be  hard  to  find 
at  300-mb.    Since  data  are  sometimes  taken  up  to  250-mb,  some  of  the 
higher  stability  values  may  be  caused  by  the  tropopause  instead  of  by 
frontal  layers.    Thus  the  value  of  this  particular  system  at  300-mb  is 
questionable.    For  example,  analysis  of  the  soundings  for  stations  906, 
909,  and  913  show  the  tropopause  between  250     and  300-mb.    This  stable 
zone  gives  high  values  of  SS  and  a  false  indication  of  a  frontal  zone. 
850-mb  1200Z  20  August  1965.     (Figure  20) 

The  major  ridge  along  the  eastern  seaboard  shows  time  continuity  in 
that  it  has  progressed  at  15  knots  over  the  eastern  states  and  25  knots  over 
the  southern  states  .    The  ridge  remains  very  close  to  the  NMC  front  through 
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the  area,  with  a  maximum  separation  of  100  miles  near  Cape  Hatteras. 
The  GG9  ridge  also  correlates  well,  except  over  Newfoundland  and  a 
questionable  area  over  the  southeastern  states.    In  this  area,  GG9  turns 
north  as  the  baroclinic  zone  separates  into  two  branches.    The  SS  ridge 
continues  across  northern  Texas,  then  northwest  into  the  Pacific,  anchored 
by  a  series  of  very  high  readings  across  the  western  United  States.    Another 
secondary  ridge,  also  in  a  baroclinic  zone,  extends  from  Delaware  west  to 
South  Dakota.    This  matches  well  with  GG9 ,  but  has  only  a  small  NMC 
counterpart  over  South  Dakota  and  has  no  time  continuity  with  any  ridge 
on  the  00Z  chart.    This  ridge  could  be  the  result  of  a  subsidence  inversion 
in  the  cool  air  behind  the  front.    The  ridge  over  Alaska  and  northwest 
Canada  is  essentially  stationary  from  00Z,  and  still  lies  in  a  baroclinic 
zone.    This  ridge  is  more  extensive  than  at  00Z,  however.    A  new  ridge  over 
extreme  northern  Canada  appears  very  strong ,  but  there  is  no  baroclinic 
zone  in  the  area . 

An  interesting  feature  of  this  chart  is  a  large  mass  of  air  with  low 
stability  which  dominates  central  Canada.    This  mass  will  be  found  to 
persist  at  all  levels,  and  seems  to  be  a  feature  of  the  atmosphere  behind 
an  advancing  cold  front. 
700-mb  1200Z  20  August  1965.     (Figure  21) 

The  ridge  line  along  the  eastern  seaboard  compares  very  closely  with 

GG9  over  its  entire  path  from  southeast  Canada  to  Montana  ,  but  west  of 

Delaware  it  is  closely  related  to  the  ridge  at  850-mb  which  was  considered 

the  secondary  ridge.    The  ridge  across  the  southern  states  which  was  strong 

at  850-mb  is  still  strong  at  700-mb,  but  is  not  continuous  across  North 
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Carolina  and  Virginia.    Considering  the  high  reading  at  station  317,  a 
higher  reading  at  station  402  would  have  allowed  a  ridge  to  intersect  the 
strong  ridge  near  station  506,  an  event  which  would  have  given  out- 
standing correlation  with  the  SS  ridge  as  shown  at  850-mb.    The  ridge  over 
extreme  northern  Canada  has  vertical  continuity  with  the  850-mb  ridge,  but 
still  does  not  lie  in  a  baroclinic  zone.    The  ridge  over  northwest  Canada  and 
Alaska  looks  very  realistic  over  its  western  portion,  but  the  eastern  half  is 
questionable  because  of  a  lack  of  data. 
500-mb  1200Z  20  August  1965.     (Figure  22) 

The  ridge  across  the  northern  United  States  is  well  defined  and  has 
consistency  with  its  GG0  counterpart.    It  also  remains  imbedded  in  the 
baroclinic  zone.    The  ridge  over  the  southeast  is  not  in  a  baroclinic  zone, 
and  appears  to  be  non-frontal.    The  ridge  over  the  southwest  desert  is  in  a 
baroclinic  zone  and  compares  well  with  GG6  in  the  area.    The  ridge  over 
Alaska  and  northwest  Canada  is  still  stationary  over  its  western  portion, 
but  enters  a  strong  north-south  baroclinic  zone  over  western  Canada,  and 
lies  very  close  to  the  GG9  ridge  in  the  area.    Another  short  ridge  over  central 
Canada  is  based  on  two  readings  and  appears  non-frontal. 

Again  at  500-mb  a  large  area  of  low  stability  exists  over  eastern 
Canada  behind  the  ridge. 
300-mb  1200Z  20  August  1965.     (Figure  23) 

The  chart  is  characterized  by  small  SS  over  the  entire  United  States 

with  the  exception  of  two  readings  of  .15  at  stations  583  and  576  which 

define  a  short  ridge  over  the  Nevada,  Idaho,  Wyoming  area.    Higher  values 

exist  over  Canada,  with  a  ridge  stretching  from  the  Great  Lakes  north,  but 
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the  tropopause  is  responsible  for  high  values  of  SS  at  stations  109,  722, 
and  943.    There  is  little  data  to  support  placing  this  ridge  entirely  inside 
the  baroclinic  zone  outlined  by  GG9,  but  contours  of  SS  are  generally 
parallel  to  the  ridge  and  trough  lines  of  GG0.    As  indicated  in  the  dis- 
cussion of  SS  at  300-mb  for  00Z,  there  appears  little  resemblance  between 
the  pattern  at  this  level  and  those  below. 
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APPENDIX  IV 
SS  Analysis  by  Method  C 
Comparison  of  SS  with  the  FNWF  Static  Stability   (<T*)  Analysis 
Technique  and  Criteria 

FNWF's  hemispheric  printout  of  (T  for  00Z  20  August  1965  displays 
areas  of  positive  and  negative  <T  (relative  to  standard  reference  values)/ 6/. 
In  some  cases  these  areas  are  large,  with  a  noticeable  gradient,  while  in 
other  cases  small  pockets  of  positive  or  negative  values  show  up. 

The  printouts  for  the  lower  four  layers  were  analyzed  for  ridge  lines 
of  0"  .    This  analysis  was  somewhat  subjective  since  no  specific  guide 
lines  were  available.    Whenever  possible,  positive  areas  were  connected 
by  a  ridge  line,  even  though  a  small  area  of  negative  values  separates  the 
positive  areas.    When  a  small  positive  pocket  was  completely  surrounded 
by  negative  values,  no  ridge  was  drawn.    The  ridge  lines  themselves  were 
smoothed  (i.e.  ,  avoiding  kinks  and  sharp  corners),  whenever  feasible,  in 
order  to  present  a  more  reasonable  analysis  .    No  distinction  was  made 
between  a  narrow,  well  defined,  positive  zone  and  a  wide,  poorly  defined 
one.    No  difference  was  noted  between  low  positive  and  high  positive  values 
either,  as  long  as  they  represented  relative  maxima. 

These  ridge  lines  were  transferred  to  Figures  15-18  and  compared  with 
SS  ridges . 
850   mb  00Z  20  August  1965.     (Figure  15) 

SS  at  850-mb  is  compared  with  the  FNWF  0"  for  the  layer  1000-775-mb. 
FNWF's    cr    is  a  meandering  ridge  line,  which  moves  in  and  out  of  the  baro- 
clinic  zone.    If  one  assumes  that  the  GG0  ridge  and  the  NMC  front  are 
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properly  located,  then   0"    is  a  rather  poor  representation  of  a  frontal  zone, 
especially  over  the  middle  west.    The  SS  ridge,  on  the  other  hand,  is 
smoother,  comparing  closely  with  GG0  and  NMC  and  lying  almost  entirely 
within  the  baroclinic  zone.    SS  and    <T  have  a  close  comparison  only  over 
the  middle  Atlantic  states.    One  factor  tending  to  produce  these  variations 
is  that   (T   considers  the  atmospheric  structure  near  the  surface  —  that  is, 
near  1000-mb.    Thus,  the  static  stability  picture  may  be  confused  by  layers 
of  air  near  the  surface  which  represent  local  terrain  features.    For  example, 
the  wide  swing  of  <T  out  over  the  Atlantic  Ocean  may  have  been  caused  by 
relatively  stable  air  over  the  ocean  east  of  Boston  —  a  not  uncommon  weather 
feature  in  summer.    On  the  other  hand,  SS  is  computed  for  the  layers  around 
850-mb,  and  this  altitude  may  be  above  low  level  inversions.    The  other 
two  or    ridges  over  the  western  United  States  and  northern  Canada  compare 
poorly  with  GG0  and  SS . 
700-mb  00Z  20  August  1965.     (Figure  16) 

Considering  the  ridges  over  the  eastern  United  States,  SS  and   C 
parallel  each  other  with  SS  generally  120  miles  on  the  warm  air  side  of  G*     . 

The  C  ridge  over  central  and  western  Canada  is  in  an  area  of  low 
SS  and  does  not  compare  well,  except  over  northwestern  Canada  and  Alaska. 
The  two  small   C  ridges  over  the  southeast  United  States  and  eastern 
Canada  are  not  verified  with  GG0  or  SS  ridges. 
500-mb  00Z  20  August  1965.     (Figure  17) 

At  this  level  both  ridges  have  less  correlation  with  the  GG0  ridge 
and  the  baroclinic  zones  than  in  lower  layers,  but  C    is  particularly  poor, 
with  only  about  25%  lying  in  a  baroclinic  zone  and  almost  no    C     ridge  to 
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compare  with  GG9.    There  are  at  least  two  locations  on  the  chart  where 
ridges  cross  GG0  ridges  at  right  angles.      CT    does  have  a  slightly  better 
positive  correlation  with  SS ,  possibly  because  the  two  quantities  have 
some  similarities. 
300-mb  00Z  20  August  1965.     (Figure  18) 

Neither    0"  nor  SS  ridges  can  justifiably  be  called  good  representations 
of  frontal  phenomena  at  this  level,  on  the  basis  of  this  chart.    As  stated 
elsehwere  in  this  project,  neither  ridge  seems  to  give  any  meaningful 
information  at  300-mb. 
850-mb  1200Z  20  August  1965.     (Figure  20) 

The   (X    ridge  over  the  northeastern  United  States  parallels  SS  and  GG9 
ridges,  generally  about  100  miles  on  the  warm  side  of  both.    However,  over 
the  middle  western  states  the    C     ridge  is  weak,  and  shows  poor  comparison 
with  extensive  GG0  and  SS  ridges  in  the  same  general  area. 

The  Cf  ridges  over  the  far  western  states  and  western  Canada  cross 
the  SS  ridges  at  right  angles  in  at  least  two  areas,  and  there  is  little  basis 
for  comparison.    The  same  is  true  of  the  ridge  over  northern  Canada  when 
compared  with  GG9  and  SS  ridges  over  Alaska. 
700-mb  1200Z  20  August  1965.     (Figure  21) 

The  highest  correlation  between     <S    and  SS  exists  at  this  level. 
This  is  a  repeat  of  the  situation  which  existed  at  00Z.    The   C    ridge 
meanders  across  the  SS  ridge,  while  crossing  the  United  States,  but 
generally  remains  within  100  miles  of  the  SS  ridge. 
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The      C     ridge  over  southwestern  Canada  is  relatively  weak  and  has 
no  correlation  with  a  SS  ridge  except  over  Alaska. 
500-mb  1200Z  20  August  1965.     (Figure  22) 

At  this  level     <T    and  SS  ridges  compare  well  over  the  northeastern 
United  States,  but  the  only  other  areas  where  comparison  is  possible  are 
over  Alaska,  a  small  area  over  California,  and  over  the  central  United 
States.    In  each  case,  the    G~  ridge  leads  the  SS  ridge,  about  100  miles  on 
the  warm  side  of  the  SS  ridge. 
300-mb  1200Z  20  August  1965.     (Figure  23) 

The  major  0*  ridge  at  this  level  compares  fairly  well  with  the  SS  ridge, 
and  very  well  with  the  GG0  ridge  over  Canada. 

The  analysis  over  the  western  states  has  no  SS  ridges  with  which  to 
compare  the    C    ridges.    The  two     <T     ridges  are  weak,  however,  and 
their  existence  is  somewhat  subjective. 

Appendix  IV  pointed  out  the  dominance  of  low  values  of  SS  at  300-mb, 
and  the    (T  printout  is  also  dominated  by  low  values ,  the  majority  of  which 
are  negative.    Both  of  these  findings  agree  with  the  conclusions  of  Gates, 
who  found  lower  values  of  static  stability  above  500-mb  /3_/. 
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